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Abstract 

We consider the production of 7*, W and Z vector bosons in hadron-hadron 
collisions in perturbative QCD. We present results from a new numerical program 
which gives a full description of the production of the vector bosons and of their de- 
cay products. At small q T the calculation includes resummation of large logarithms 
and non-perturbative effects. The resummation is matched with the full 0(as) 
calculation. In addition, the program correctly reproduces the known 0(ots) cross 
section when integrated over q T . Besides presenting results for W and Z production 
at the Tevatron, we also review constraints on the non-perturbative functions using 
fixed target data on lepton pair production, and make several observations on this 
topic. 



1 Introduction 



The aim of this paper is to give a complete description of the leptons coming from the decay 
of a vector boson produced in hadron-hadron collisions. In addition to the intrinsic interest 
of obtaining a complete description of vector boson production, a precise measurement of 
the W mass at a hadronic collider will depend on accurate theoretical information about 
its production properties. Since it is the leptonic decay products of the vector bosons 
which are actually observed, it is important to include the decay so that experimental 
cuts can be implemented. Leptons produced at large transverse momenta can come from 
vector bosons having either large or small q T , where q T is the transverse momentum of 
the vector boson. It is therefore essential to use a formalism which is valid for both large 
and small transverse momentum of the vector bosons. 

The production of vector bosons at large and small q T has been extensively consid- 
ered before. Vector bosons at large q T attracted early interest because they provided 
discrimination between the naive parton model (which predicts a limitation on q T ) and 
an underlying field theory (which predicts events at large q T ). The transverse momentum 
of a vector boson recoiling against one parton was calculated in Ref. [|l|. The 0(a 2 s ) 
calculation of vector boson production at large q T was initiated in Ref. [Q] and completed 
in Refs. and [|J . These papers on the q T distribution give no information on the distri- 
bution of the leptons into which the vector bosons decayed. This deficiency was remedied 
in Ref. H for the 0(as) case and in Refs. || [7j for the 0(a|) case. 

However the bulk of the data is not at large q T . At small q T order by order in pertur- 
bation theory we encounter large logarithms of Q 2 /q 2 , where Q is the vector boson mass. 
For example, at small q T the leading term in the cross section is of the form 

^~^ln^. (1) 

dq 2 q 2 q 2 

The logarithms can be resummed to give a Sudakov form factor || ||. The resummation 
has a simple exponential form after transformation to the impact parameter, b, which is 
the Fourier conjugate of q T 0. The necessary coefficients for the inclusion of higher order 
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terms in the resummation were calculated in Ref. |10[ using the results of Ref. H. The 



formalism which we shall use was written down in Ref. JTTJ using techniques developed 



earlier for back-to-back jets [[L2[]. Numerical results of resummed calculations have been 
reported in Refs. [0, [L4j]. 

In order to have a complete description of the transverse momentum one must match 
the theoretical results at large and small q T . The matching of the vector boson cross 



section including 0(ag) has been considered in Ref. The matching including the full 



0(a.g) calculation appropriate at large q T is given in Ref. |T5|]. The papers in Ref. |T6 
have combined fixed order calculations with a parton shower approach. 

One of the most interesting results of the resummation procedure is that for large 
enough vector boson mass perturbation theory is valid even at q T = 0. In fact, using the 
saddle point method it is found || |TT| that the Fourier transform integral to recover the 
transverse momentum distribution is determined (at q T = 0) by an impact parameter of 
order 

where « = 16/(49 - 2n f ). Choosing Q = M w = 80.33 GeV, A = 250 MeV and the 
number of active flavours n/ = 5 we find that the saddle point is at the position 

6gp « 10A . (3) 

Thus for W and Z production the integral is dominated by values of b which are at the 
borderline between the perturbative and non-perturbative region. Detailed predictions 
therefore depend both on the perturbative Sudakov form factor and on a parametrization 
of the non-perturbative part of the form factor, to be extracted from data. The effect 
of non-perturbative terms on the vector boson q T distribution has been considered in 
Refs. 00,0- 

A treatment of the production of vector bosons including both the resummation at 



small transverse momentum and the decay kinematics has been given in [ I9fl . In this 



paper we carry the analysis further and correct minor mistakes in Ref. [[Tj|. A more 
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complete theoretical description could be obtained if the order a| were included at large 
q T , (i.e. by extending the results of Arnold and Kauffman [IS| to include the decay of the 
vector boson). 

In Section |2| we give a brief summary of the formalism which we use for the theoretical 
description of vector boson cross production. Section [| presents results for the 0(as) 
g T -integrated cross-section. Section ^ contains our numerical results. Besides discussing 
W and Z production, we also review the available information on the non-perturbative 
parameters which are important at low q T , and point out several issues which have not 
been addressed previously in the literature. Our g T -dependent predictions are obtained 
and verified using two independent programs based on the formalism of Section |2|. For the 
g T -integrated distributions we have implemented a separate program based on the results 
of Section || Our conclusions are presented in Section |[ while details of the formulae 
used in our numerical programs are given in two appendices. 



2 Resummed cross section at hadron level 

In the Collins-Soper frameQ |2(J the general expression for the resummed differential cross 
section at hadron level may be written in the form 

da{AB -> K(-> ll')X) 1 Q 2 



dq 2 dQ 2 dy d cos 9 # 2 s NnS (Q 2 - M 2 ) 2 + M 2 T 2 V 



x 



Y r (q 2 T ,Q 2 ,y,9) + Y f (q 2 T ,Q 2 ,y,9, 



(4) 



In the above, N = 3 is the number of colors, v^S* is the total hadron-hadron center-of-mass 

energy, while 9 and <fi refer to the lepton polar and azimuthal angles. The functions Y r 

and Yf stand for the resummed and finite^ parts of the cross section, respectively. They 

are defined in the following subsections. 

1 The Collins-Soper frame (defined in the Appendix |a|) is the rest frame of the vector boson with a 
specific choice for orientation. 

2 The "finite" part is integrable as — > and contains no distributions. 
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2.1 Resummed part 



The resummed part of the cross section is given as the Fourier integral over the impact 
parameter b, 



Y r (q 2 T ,Q 2 ,y,6) = -L [bdbJ (q T b)Y,'Fa b P (Qib,^x B ) 

(1,0 

x Wab(Q, K, 0)f' a/A (x A , fi(b*))f{, /B (x B , 11(h)) 



(5) 



In this expression the function F NP represents the non-perturbative part of the form 
factor. Its specific form, as well as the definition of the variable 6* will be described 
below. The variables x A and x B are given in terms of the vector boson mass Q and 
rapidity y as 



' \ -^Lexp(y) , x B = -Q=exp(-y) . (6) 



S 



Note that 



t — x A x B — Q I ' S 



(7) 



The modified parton structure functions /' are related to the MS structure functions 
/ by a convolution, 

- 1 dz 



fa/A( x A, H) = f — C ac( — , V)fc/A(z, V) 
„ Jxa Z Z 



(8) 



where (a, b 7^ g) [p 



^ag{Z,H) — — J-R 

111 



2tt 
2z{\ - z 



,7T 



(9) 
(10) 



In the above the colour factors are Cp = 4/3 and Tr = 1/2, while the prime on the sum 
in Eq. (|5]) indicates that gluons are excluded from the summation. 

The function W from Eq. (JD can be written in terms of the the Sudakov form factor 
S as 

W ab (Q,b,8) =exp[S(b,Q)]Hj°\9) . (11) 



5 



The exact definition of the Sudakov form factor will be discussed below. The function 
H(°\ which includes the angular dependence of the lowest order cross section and coupling 
factors, is defined in Appendix [A|. 



2.1.1 Sudakov form factor, large and small b 

In the formalism of Ref. [|11|] the Sudakov form factor S(b, Q) is given by 

rQ 2 djl 2 



S(b } Q) 



b 2 /b* /i 2 



(12) 



with bo = 2exp(— 7^) 1.1229. The coefficients A and B are perturbation series in a,s, 



2tt/ 



(13) 



where the first two coefficients in the expansion are known [ 10[ : 
A {1) = 2C F , 

A« = 2 C F (iV(|-^)-H TRlv 
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B {i) 



-3C f 



5 (2) = c^-l-m^+c.N^-^+eas) 



4 

+ C F T R n f (-j - -it 2 



(14) 



In the above we take rif to be the number of quark flavors active at the the scale at which 



0£g is evaluated. In addition, the formalism of Ref. [JTTJ] requires that the scale at which 
the parton distributions are evaluated in Eq. (Bp is 



fj,(b) = b /b . 



(15) 



Equations (|I2"D through QTSJ) should be compared with the exact first order results JTS 
for the Sudakov form factor, 



and for the scale (i(b), 

fl (b) = Qexp{-J o Q ^[l- J (bx)]} . (17) 

Formally, the integration over b in Eq. (||) is from to oo. However, as b approaches 
1/A, the coupling as becomes large and the perturbative calculation of the form factor 
S is no longer reliable. This region is effectively removed from the integral by evaluating 
W and the parton structure functions at 

(18) 



lim) 



which never exceeds the cut-off value b \ im . The large b part of the Sudakov form factor is 
provided by the function F NP which parametrizes the non-perturbative effects [JTOL |TT[ . 
The specific form of F NP , as well as our results with particular choices for non-perturbative 
parameters, will be discussed in Section [|. 



The formalism of Ref. |TTJ leaves open the question of small b. The small 6-region does 
not contribute large logarithms, but a correct treatment is important to recover the total 
cross section after integration over q T . The lowest order expression of Eq. ( Pp had the 



property that S — > as b — > 0, which is lost in Eq. fll2f) . In Ref. |2]J it has been suggested 
that one make the replacement in Eq. (|T!|) 

{ b } { b } [1 + 6§/(6 2 Q 2 )] ' { } 

This ensures that the scale never exceeds Q by adding power suppressed terms of order 
(Qt/ Q) 2 - 

However, we shall use a slightly more sophisticated treatment which ensures that the 
first-order fixed-coupling result is correctly reproduced. We define scales A and \x by 

/ — In ^ = / — In ^ 1-Jo te , (20) 

J\ X X Jo X X 

Qi ± = j Qd ± {l - Uix)] , (21) 
, x Jo X 
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which results in 



Kb) 
A(6) 



Qexp{-| o Q ^[l- J (bx)]} , 



Qexp{-[2 / Q ^ln(^)[l- J (6x)]l 5 } 



(22) 
(23) 



10 X X 

Figure shows A and /J plotted for Q = 5 and 100 GeV. At large b ^> 1/Q where the 
resummation is mandatory, we have that 

6 



A, a* 



(24) 



in accordance with Eq. ( |l~2"l) and the procedure of Ref. |Tl|- In addition we have that 
n(b), A(6) < Q, where the equality is true for b = 0. 

Therefore, instead of Eq. (|12"D for the resummed form factor we shall use 

dfi\_ (Q 2 \ A ,_ ^ , f Q2 dfi 2 



5(6, Q) 



(25) 



For fixed coupling constant this expression is exactly in agreement with the lowest order 
result of Eq. ([16]) , but also preserves the good features of Eq. fll2"|). The exponential of 
the Sudakov form factor with this prescription is shown in Figure |2] for Q = 5, 10 and 
100 GeV. 



2.2 Finite part 

The finite part of the cross section Eq. (|j) is defined as 

71 V Jx A Z A Jx B Z B ab . =1 

x /^(Q 2 , z A , * B , 0, 0) / oM (-, 0)/ fc /fl(-, Q) • (26) 

To order 0{a$) we only need the function w- l \ which is given by the difference of the 
parts derived from one parton emission, and the pieces which have been removed from the 



a s (Q) 
2tt 
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cross section in the factorization or resummation procedure (H^ and respectively). 
Therefore, we can write 

Ra b \Q 2 ,z A ,z B ,9,(f)) = Q 2 Hi 1 b \Q 2 ,z A ,z B ,9,4 > )6(s + t + u-Q 2 ) 

- ^ b \Q 2 ,z A ,z B ,6) . (27) 

can further be separated into parts which are divergent (S^) or integrable as q 2 — > 
i.e. 

H2? = H2? + g$. (28) 

S'W can then be combined with into a function £W, so that the whole residue 
now in the form 

R ( £(Q 2 ,z A ,z B ,e,<f ) ) = Q 2 H^ b > (Q 2 ,z A ,z B ,e, ( f ) ) 5(s + t + u~Q 2 ) 

+ ££(Q 2 ,* A ,* B ,0) , (29) 

has the property of being integrable as q 2 goes to zero. Explicit expressions for functions 
H^ 1 ' and £W, as well as the relation of Q 2 , s, t and u to z A and z B , are given in Appendix 



3 Integrated cross section 

With our definitions, after integration over g T ,cos# and 0, and dropping 0{a 2 s ) terms, 
we recover exactly the order 0{a,s) cross section, which can be written in the form 

da(AB — > VX) 1 Q 2 16 



dQ 1 dy 2 8 NttS (Q 2 - M 2 ) 2 + M 2 T 2 V 3 

&s(Q) f 1 dz A f 1 dZj: 



^ ' f'a/A( X ^i Q)fb/B( X A, Q) 



o,6 



+ — XI ^(^a, z B )f a / A {—, Q)h/B(—, Q) 

ZTT Jx A Z A Jx B Z B „ , Z A Z R 



.(30) 



/x A Z A Jx B Z B ab Z A 

In this case the modified parton distribution functions /' are defined as (a, b ^ g) 

L^A./i) = Yl [ — D ac( — ,^)fc/A(z^), (31) 
n JXi. Z Z 
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with coefficients 



2tt 



1 — z 111 



1 - z 



D ag (z,n) 



+ {1 + Z 

2tt 



/ln(l 



,7T 



V 1 - z 



+ ( __ 4 ) 6(1 -z) 



T, 



R 



2z(l-z) + \z 2 + (l-z 



2(1 -z) 



+ z 



(32) 
(33) 



As before, the prime on the sum indicates that gluons are excluded from the summation. 
The functions X ab are given as (X qq = X qq , X qg = X qg , X gq = X gq ) 



X qq {z A , z B ) 

-Xqg(z A , Z B ) 

X gq \z A , z B ) 
where we defined 



Gq ' {ZA ' Zb) ((1 - z A )(l - Z B 



H~ F qq (z A , z B j , 



Gqg (z A , Z B ~) 

G qg (z B , z A ) 



1 



1 — Zk 



1 — Z r 



Gq q {z A , Z B ) 



c F 
c F 



2(1 + z A z B )(z 2 A + zl) 

(z A + z B ) 2 
2(1 + z A z B )(zj + zl) 

(1 + Z A )(1 + Z B ) 

2z A (l + z A z B ) 
(1 + z A )(z A + Z B ) 



' hZAZB)2zl - + (l-z A z B ) 2 



z A + Z B 



(34) 
(35) 
(36) 

(37) 
(38) 
(39) 



The definitions of the single and double "plus" distributions used in the above are given in 
Appendix |A]. We also remind the reader that results in this section are expressed in terms 
of MS scheme structure functions. Similar expressions in the DIS scheme are presented in 
Ref. Using the techniques of Ref. |22 we can further integrate this to the standard 

result for the total cross section []23| in the MS scheme. 

3 We believe that the qg contribution in Ref. |22] is in error. 
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4 Results 



As already noted in the introduction, the motivation for this work is that a precise mea- 
surement of the W mass at a hadronic collider will depend on accurate theoretical in- 
formation about its production properties. However, before presenting our results for W 
and Z production, we would like to address the following issues which in our opinion 
have not been adequately discussed in the literature: the determination of the form of 
the non-perturbative function from the low-energy Drell-Yan data, the dependence of the 
results on the choice of bn m and the matching between low and high q T . 



4.1 Determination of F 

The unknown function F NP from Eq. ([|) has a general form [1 I 



F ti (Q,b,x A ,x B ) = exp<- 



h Q {b) In (7^-) + hi(b, x A ) + hj(b, x E 



(40) 



where the functions h are not calculable in perturbation theory and therefore must be 
extracted from experiment. On general grounds, we expect that hi — > as b — > 0, so 
that the g T -integrated cross section is unchanged. On the other hand, the parameter Q 
is completely arbitrary. 

The first attempt to obtain F NP from experiment was made by Davies et al. (DWS) 



in Ref. |L4|. There the functions h were approximated by 

hq{b) = g 2 b 2 , 
hi(b,x A ) +hj(b,x B ) = gib 2 , 



(41) 



since the observed q T distribution at low Q was approximately gaussian in shape. Using 
the Duke and Owens parton distribution functions pi |, DWS determined parameters g\ 
and g 2 from E288 |25| data with y/S = 27 A GeV, and also from R209 |26| data. The 
resulting values (with a particular choice of Qq = 2 GeV) were 



g x = 0.15 GeV 2 , g 2 = 0.40 GeV 2 



(42) 
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The cut-off vafue of b u m from Eq. (|i~8| ) was chosen to be 0.5 GeV -1 . This parameter set 
yielded a good agreement of theory with R209, as well as with E288 data for Q < 9 GeV 
mass bins. However, for Q > 11 GeV theoretical expectations were unacceptably far 
above the data. 

Motivated by the fact that the production of vector bosons at Fermilab Tevatron 
(\/S = 1.8 TeV) involves values of r = x A x B which are significantly lower than those 
considered in Ref. [rj|, Ladinsky and Yuan (LY) JT7] reinvestigated the form of the non- 
perturbative functions from Eq. (|40|) . Using CTEQ2M parton distribution functions [pTfl , 
they showed that the DWS form of non-perturbative function (Eqs. (|41] and (fl~2^ )) no 
longer agrees with R209 data for 5 GeV < Q < 8 GeV. In order to improve theoretical 
predictions, LY postulated the r dependence for functions hi, so that 

M fo ) = 92b 2 , 

hi(b,x A ) + hj(b,x B ) = g x b 2 + g x g 3 b In (— ) , (43) 

where r is arbitrary parameter.^ Choosing r = 0.01, Qo = 1.6 GeV, and b ij m = 
0.5 GeV -1 , these authors determined the non-perturbative parameters by comparison 
of theory to R209 data |2BJ in the range 5 GeV < Q < 8 GeV, to ^/S = 27 A GeV E288 



data |25[ in the range 6 GeV < Q < 8 GeV, and also to CDF Z data p8fl . In particular, 
LY explicitly showed that values of 

9l = O.Ut°of 3 GeV 2 , g 2 = 0.58^ GeV 2 , ^3 = -1.5lg;l GeV" 1 , (44) 

provide a good agreement of theory with CDF Z and with R209 data (for 5 GeV < Q < 

8 GeV), as well as with CDF W data [29]. They furthermore noted that the parameters 

of Eq. fl44|) give results which are in agreement with E288 data and with R209 data for 

11 GeV < Q < 25 GeV. 

The problem with using the Drell-Yan data from the fixed-target experiments is that 

the overall normalization of the cross section is uncertain. For example, the E288 data 

4 We note here that, unlike a simple gaussian, the particular functional form of the non-perturbative 
function in Eq. ( |43| ) is not always positive. 



12 



has a stated normalization uncertainty of 25% [pSfl . Since smearing function F simply 
shifts the q T distribution between the low and high q T regions, and since the bulk of the 
data is in the low q T region, it is clear that using data with wrong normalization will 
affect the non-perturbative parameters. Therefore, it is necessary to establish correct 
normalization for fixed-target experiments before trying to determine the shape of F NP 
from their data. 

A consistent way of determining overall normalization is to compare theoretical pre- 
dictions for g T -integrated cross section (to a given order in as) with experimental results. 
We illustrate that procedure for E288 |25j and E605 |3(J experiments.^ 

For E288 (y/S = 27 A GeV) we used S da 2 /dy/rdy distributions for the y/r bins with 
Q < 9 GeV.Q In order to achieve agreement with experiment (see Figure |3|), we found that 
theoretical results had to be rescaled down by a factor K = 0.83 ± 0.03.0 On the other 
hand, for the invariant cross section versus q T data in the range 5 GeV < Q < 9 GeV and 
q T < 2 GeV, the central values of LY parameters yield the best x 2 °f about 7.7 /do f for 
K = 0.75, while ^-factors of 0.80 and 0.83 lead to x 2 of about 11.9/ do f and 18.2 / do f, 
respectively. These results, shown in Figure |J indicate that parameters of Eq. ([14]) 
overestimate the E288 transverse momentum data. We are therefore unconvinced that 
the parameters of Eq. (|44]) give the best possible fit to the data. 

For E605 (y/S = 38.7 GeV) we used all available y/r bins for S do 2 / d^frdy distribu- 
tions, and determined .fT-factor of K = 0.88 ± 0.02.Q The agreement between the theory 
and the data is illustrated in Figure |^. The corresponding invariant cross section distribu- 
tions obtained with central values of LY parameters are shown in Figures || and [f[ From 

5 Unless otherwise stated, all results described in this section are obtained using CTEQ2M parton 
distribution functions, which facilitates comparison to previous work p7| . 

6 The y/r bins with Q > 11 GeV were discarded because of the low statistics. 

7 We define the if-factor as K = experiment / theory . We have chosen to change the normalization 
of the theory, rather than shifting the experimental data, despite the fact that it is the data which is 
subject to a normalization uncertainty. 

8 E605 experiment has a stated normalization uncertainty of 15%. 



13 



Figure |6| it can be seen that theoretical predictions are in good agreement with experimen- 
tal results for the bins with Q < 9 GeV. However, for bins with Q > 10.5 GeV (Figure 
|7|) the theoretical distributions do not match the distributions obtained by experiment. 

The LY functional form of F NP given in Eq. (|43|) implies that theory should be able 
to describe data with different values of r. Our results certainly do not support that 
statement, and we believe that the form of the non-perturbative function remains an 
open question. 



4.2 Choice of &r im 

Another issue which has not been addressed in the literature is the dependence of results on 



the choice of the cut-off value bu m from Eq. ([T^) . In the original work of Ref . |14]] b u m was 
taken to be 0.5 GeV -1 because the structure functions were not defined for scales less than 
2 GeV. However, this choice is arbitrary. Any change in b n m (in a reasonable range around 
0.5 GeV^ 1 ) should be compensated by a change in non-perturbative parameters describing 
F NP , so that an equally good description of experimental data is always achieved. 

In order to verify that statement, we attempt to reproduce several different sets of 
transverse momentum distributions, with b\ im chosen in the range from 0.3 GeV" 1 to 
0.7 GeV~\ which one may consider to be reasonable.^ We have chosen ^-distributions 
that include R209 dajdq% data (5 GeV < Q < 8 GeV, with q T < 3 GeV), E288 Eda/d 3 p 
data for Q = 5.5, 6.5, 7.5 and 8.5 GeV (with q T < 2 GeV and i^-factor of 0.83), and also 
E605 Eda/d 3 p data for Q = 7.5,8.5, 11 and 12.5 GeV (with q T < 2 GeV and ^-factor 
of 0.88). Since these data were obtained within a narrow mass range we can neglect the 
Q-dependence of Eq. fl40D, and adopt a simple gaussian form for F NP , 

F NP = exp (-g b 2 ) , (45) 

where g is an effective parameter, different for each of the above q T distributions. In this 



9 Since CTEQ2M structure functions are denned for scales greater than 1.6 GeV, and since fi(b) ~ bo/b 
for large b, the upper limit of 0.7 GeV -1 is still a possible choice for b n m . 
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way, the task of finding the form of the non-perturbative function is reduced to a simple 
one-parameter fit. In Figures I8J, and ITU] we show the best x 2 /dof obtained by varying g 



with different values of bn m , for R209, E288 and E605 data sets, respectively. It is obvious 
that for R209 data (where y/r = 0.105) an equally good fit can be obtained regardless of 
the value chosen for b lim , in part because there is only a small number of data points of 
limited statistical precision. 

However, for E288 and E605 data, which involve larger values of \fr (for E288 \fr = 
0.201,0.238,0.274 and 0.311, while for E605 ypF = 0.194,0.220,0.284 and 0.323), the 
choice of b u m makes a considerable difference, which is not satisfactory from the theoretical 
point of view. These results might indicate that perhaps a pure gaussian form of the non- 
perturbative function is wrong.Q In any case, it is clear that this problem requires a 
further study. 

In Table |TJ we also note the range of values for g for which results shown in Figures 
H through [IJJ were obtained. If the gaussian form of the non-perturbative function were 
correct, then the values of g, which increase for increasing Q in the E288 and E605 sets, 
would indicate that the effective coefficient in front of In (Q/(2Qo)) in Eq. (|40D should be 
positive, and furthermore that the r dependence of F NP is not large. However, the values 
of g obtained for R209 data do not seem to support that observation. 

4.3 Matching of low and high q T regions 

The resummation formalism is expected to give a good theoretical description of vector 
boson production in the low q T region (g| <C Q 2 ). On the other hand, conventional 
perturbation theory provides a good approximation in the other regime, for q 2 ^> Q 2 . 
The necessity of matching low and high q T regions has already been discussed in Ref. |15] 



for the W and Z production, where it was shown that the proper matching of the pure 
perturbative and resummed expressions reduces theoretical errors. However, matching 



10 There is no reason in principle why the functions h from Eq. (^0|) should be limited to a quadratic 
dependence on b. 
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will never be perfect since resummation introduces higher order terms which will not be 
cancelled at large q T in any finite-order calculation of Yf from Eq. (^). For this reason, 
the conclusion of Ref. |15] was that one should prefer the ordinary perturbation theory 
result once the resummed part Y r becomes negative. For the W and Z production with 
= 1.8 TeV this happens at q T ~ 50 GeV. In Ref. |15] it was also shown that matching 
works well with the second-order calculation of Yf. For example, at q T of about 50 GeV 
the mismatch between the q T distribution calculated using resummation plus the 0(a 2 s ) 
calculation of Yf and the one calculated using conventional second-order perturbation 
theory was of the order of 10%. 



In Figure [11] we show our results for the W production, obtained with gaussian and 
LY form of F NP , compared to the 0(cts) perturbative result. At q T of 50 GeV the 
mismatch between Y r plus the 0(as) calculation of Yf and perturbation theory is about 
50%. For the 7* production at fixed target experiments the same problem (but even more 



acute) is illustrated in Figure [12] for E288 experiment (with 5 GeV < Q < 6 GeV and 
—0.27 < y < 0.33). These two figures clearly indicate the necessity for extending the 
results of Ref. fl5|] to include the decay of the vector boson, if a complete theoretical 
description is desired. Nevertheless, we emphasize that results presented in this paper 
are obtained for the low q T region, where the resummation formalism plus the 0{ots) 
calculation of Yf should be adequate. Note also that it is the region of low q T which is of 
interest for the W mass measurement. 



4.4 Results for the W and Z production 

Once the high statistics data on the vector boson q T distributions at 1.8 TeV become 

available, one can extract the effective form of F NP from the Z data, and use that 

information to obtain an accurate theoretical prediction for the various g T -dependent 

distributions of the W boson.[[]] This in turn should allow a precise measurement of the 

11 Taking current values of vector boson masses pi] , at VS = 1.8 TeV for Q = M W (M Z ) the value of 
y/r is about 0.045(0.051). 
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W mass. 

For the data available at present ^] the statistics is low, which limits the predictive 
power of the resummation formalism. In order to illustrate that, we again adopt a simple 
parametrization of F NP given in Eq. fl45|), choose b\\ m = 0.5 GeV -1 , and vary g in an 
attempt to obtain a good fit to the W = W + + W~ and Z da/dq T data.0 As shown 
in Figure |13], due to large statistical errors almost any choice of g in the range between 
of about 1.5 GeV 2 to about 5 GeV 2 yields an acceptable description of both data sets. 
We conclude that it does not make much sense to use these data for determination of 
F , and that data with much higher statistics are needed before any firm theoretical 
predictions can be made for the W and Z production. Nevertheless, we still observe 
that the effective value for g obtained from the Z data tends to be smaller than the one 
obtained from the W data. This is in disaccord with what one would expect from the 
fixed target data (see Table [I]), and may again indicate that pure gaussian form for F NP 
is not correct. It may also indicate experimental biases introduced by the selection of two 
isolated leptons in the Z sample. 



In Figures [LJ and [lg we show our results for the W and Z q T distributions. These re- 
sults are obtained with F NP given in Eq. (f45"[), with g = 3.0 GeV 2 and b \ im = 0.5 GeV -1 , 
using several different parton distribution functions. Besides illustrating the as depen- 
dence of da/dq T j^\ these two figures also show that the fit to the data is as good as the 
one obtained in Ref. jfTTj] , even though we used much simpler functional form of F . 

We now consider briefly the import of these results for the measurement of the W 

mass. Figure [H] shows the transverse mass tjit of the lepton pair, obtained with the two 

different choices of F NP . To first order the transverse mass is insensitive to the transverse 

motion of the W, and because of that the ttit distribution is largely independent of the 

non-perturbative parameters. 

12 We assumed BR(Z -> e+e") = 0.033 and BR(W+ -> e+u) = 0.111, as was done in @ and g|, 
respectively. 

13 For CTEQ2M 0, MRSR1 and MRSR2 §1, we used a s (M z ) of 0.110, 0.113 and 0.120, respectively. 
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As the luminosity of the Tevatron is increased the number of interactions per beam 
crossing will increase, leading to a degradation of the missing energy resolution. Therefore 
the measurement of the Jacobian peak in the lepton transverse momentum will become a 
competitive method of measuring the W mass. Figure |j] shows the expected transverse 
momentum distribution of the electron from W~ decay. The width of this distribution is 
broader than the transverse mass distribution and the dependence on the non-perturbative 
functions is larger. A quantitative estimate of the size of this dependence will have to 
await a reliable extraction of the non-perturbative parameters. 



5 Conclusions 

In view of the large number of W and Z bosons to be expected in Run II at the Tevatron 
we have returned to consider their production and decay in hadronic collisions. We 
have provided a description of vector boson production which not only gives a correct 
description at small q T , but also reproduces the correct formula for the g T -integrated cross 
section. In addition we have included the decay of the vector bosons so that experimental 
cuts can be included. 

In the course of our numerical work we have raised several issues which have not been 
adequately addressed in the literature. The analysis of the low energy experiments needs 
to be repeated, using the g T -integrated data to fix the overall normalization, before any 
attempt to determine the form of the non-perturbative function is made. Furthermore, 
this analysis should include all low energy experiments for which g T -dependent distribu- 
tions are available. On the other hand, for the W mass measurement the effective form 
of the non-perturbative function can be extracted from the Z data. This would eliminate 
uncertainties related to the determination of F NP from the low energy experiments. 

Besides the form of the non-perturbative function, there are also other theoretical 
problems which have to be resolved. In particular, we have shown the necessity for 



extending the results of Ref. |0 to include the vector boson decay, if a more complete 
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theoretical description of the leptons coming from that decay is desired. As shown in 



part of Eq. (|j) should yield satisfactory results for matching of low and high transverse 
momentum regions. However, it is not quite clear whether such a calculation would 
entirely solve the problem of getting a complete description of the q T distribution for the 
7* production in the low energy experiments. 



This work was supported in part by the U.S. Department of Energy under Contracts No. 
DE-AC02-76CH03000. R.K.E. and D.A.R. thank the CERN theory group for hospitality. 

A Details of the formulae 
A.l Couplings 

In this appendix we document the results for functions H^°\ and which appear 
in Eqs. ([11]) and (|29|). These functions can all be separated into parts which are even and 
odd under parity, e.g. we can write 

H (o) = H (o)+ + H (o)-^ Consequently, we first define 
"plus" and "minus" quark-quark and quark-gluon couplings as 



Ref. [r5||, for the W and Z production at Tevatron the 0(a 2 s ) 



calculation of the finite 



ACKNOWLEDGMENTS 



vt q > = \v m >\ 2 {gl + gl){fl + fl), 

V w , = \V qq/ \ 2 (g 2 L -gl)(ft-fl), 



(46) 
(47) 



and 



v g + ff = Ei^r^ 2 +^)(/ 2 +/ 2 ), 



(48) 



v qa = Y.\v< M >\ 2 i9l-9i)Ul-fl), 



(49) 
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where g L , g R , f L and / R are listed in Table |2|. The coefficients V qq > are elements of the 
Cabibbo-Kobayashi-Maskawa matrix for the W production, and are equal to 5 qq > in the 
case of Z or massive photon 7*. 

In the case of production, which can proceed through the exchange of either a Z 
or a 7*, the above expressions need to be modified by making the replacements (Q e = — 1, 
Q u = 2/3, and Q d = -1/3) 

{9l + gl)Ul + fD - (gt + gl)(fZ + fl) + ^Q}Qlx2(Q 2 ) 

+ 2e 2 Q f Q e (gJ L + g R f R + g h f n + g K f L ) X i(Q 2 ) , (50) 

{9l-gl){fl-fD - (gt-glKfl-fl) 

+ 2e 2 Q f Q e (gJ h + g n f n - gj n - gMxi(Q 2 ) , (51) 



with 



xm = , (52) 



(Q 2 - Ml) 2 + M 2 T 2 , , 

X2(Q 2 ) = y - ^ (53) 

In these expressions e 2 is the electromagnetic charge which is taken to run (down from 
its value at the Z coupling) using the one loop electromagnetic (3 function, 



da a 



din Q 2 37r L ^ 

At Q = 5 GeV we find that a ~ 1/133. 

In our numerical work we choose four input parameters 14 

G F = 1.16639 x 10" 5 GeV" 2 



14 The boson widths are fixed at their measured values 

Y w = 2.07 GeV, T z = 2.49 GeV . 
These, as well as all other parameters, are taken from the Review of Particle Properties 
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M z = 91.187 GeV, 
a(M z ) = (128.89)" 1 , 
M w = 80.33 GeV . 

In terms of these parameters we can derive the p parameter 



(55) 



1 - 



na(M z 



y/2G F M? 



w \ 



1.00654 



(56) 



which enters in the couplings (see Table i n the improved Born approximation: 



9w 
9l 



x w 



AV2M^G F , 
V2G F M 2 z p , 



1 - 



Mw 



(57) 



A. 2 Matrix elements 

We start by considering the lowest order process for the production of a vector boson of 
mass My and width Yy (the momenta are shown in brackets), 



?(pi) + qipz) -> Kfa) + 1(h) ■ 



(58) 



The matrix element squared for lowest order process averaged (summed) over the 
initial (final) spins and colours (N = 3) is given by 



Q 4 



-H, 



(o) 



AN (Q 2 - Myf + MyTy « ' 



(59) 



where 



H, 



(o) 
qq 



(pi • & 2 p 2 • &i + pi • fci p 2 • fa) + V q g (pi -k 2 p 2 -fa-p\- fa p 2 ■ k 2 



(60) 
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withF] 



It is convenient to express the above matrix element in terms of angular variables in 
the Collins-Soper (CS) frame, which is defined by 

ki = — (1, sin0sin#, cos sin 6 1 , cos 9) , 

k% — y (1, - sin sin 0, - cos sin 9, - cos 9) , 

Pf = -^^(1,0, -sin /5, cos/?) , 

= - {U ~ 2 q 2) (l,0,-sinA-cos/3) , (61) 

tan/3 = ^ . (62) 

In the terms of the vector boson transverse momentum (q T = |q T |) and rapidity (y) in 
the lab frame, we also have 

qf = (M T cosh y, q T , M T sinh y) , (63) 

where M|. = Q 2 + g| and Q 2 = g 2 . We further define several functions which determine 
the angular dependence as 

Co = 1 + cos 2 9 , 
A = ^(l-3cos 2 fl) , 
Ai = sin 20 cos , 
A% = - sin 2 6* cos 20 , 
^4 3 = 2 cos , 

A4 = sin 9 cos . (64) 
Using these angular variables Eq. ( j60|) may be written as, 



H$(9)=V+£o + V7,A 3 . (65) 



15 In lowest order the vector boson is produced at zero transverse momentum and hence = 0. 
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Note the relation between the qq and qq processes, 



(66) 



We next consider H^ 1 ' and S*- 1 -* (needed for the finite part of Eq. (f|)), which are derived 
from the matrix element squared for processes involving one parton emission, and from 
the pieces removed from the cross section in the resummation procedure, (cf. Eqs. 
and (p9f) ). It is sufficient to calculate results for the two-to-three processes 



?(Pi) + Q(P2) 

q{pi) + gfa) 



l(h) + l{k 2 ) + g (k 3 ) 
Kh) + l{k 2 ) + q '{h) 



(67) 
(68) 



as all other processes are determined by the crossing relations. The invariant variables 
for the above processes are 

s = ( Pl +p 2 ) 2 , 

g 2 = Q 2 = (k 1 + k 2 ) 2 , 

t = (pi - qf = (P2 - hf , 

U = (P2 - ?) 2 = (P - & 3 ) 2 ■ 



The matrix element squared can be put in the form 

Q 2 



9 



2N (Q 2 - My) 2 + MyT 



H. 



(i) 



2 p2 a& ' 



(69) 



(70) 



v 



with 



(i) 
qq 



8C f 



til 



V, 



qq 



+ v, 



#(1) 

99 



99 

8T, 



(pa • ^i) 2 + (pi • fe) 2 + (p 2 • A; 2 ) 2 + (Pi • fei) 
(/-A:i) 2 + (pi-A: 2 ) 2 -(p 2 -A ;2 ) 2 -(pi-A: 1 ) 2 

(A* • h) 2 + (p! • fc 2 ) 2 + (fc 3 • A; 2 ) 2 + ( Pl • feO 



(71) 



9.9 



+ v, 



9.9 



(A* • h) 2 + ^ • A; 2 ) 2 - (fc 3 • A; 2 ) 2 - ( Pl • ^) ; 



(72) 
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Again, using Eq. fl6l|) the functions can be expressed in terms of invariant variables 
and the angles in the CS frame 

H ab H$(Q 2 ,t,u,6,(j)) . (73) 

Before writing down the expressions for and E^, we clarify the relationship 
between different variables. The Mandelstam variables for the parton subprocess can be 
expressed in terms of the integration variables in Eq. ( pB[) (z A = x a /£a and z B = a; B /£ B ) 
as 

Q 2 

s = , 

t = (?7 - Za) , 

u = _^! (r? _^ B ) ? (74) 



z, 



B 



where M T = yQ 2 + q 2 and r\ = M T /Q. We further note that if Q 2 = s + t + u and 
q 2 = ut/s one can derive 

_ ^ = 1 + ^B ^ 
Q 4" ^B 



so that Q 2 {= q\j{j] 1 — 1)) can be expressed in terms of q 2 , z A and z E 

r(l) _ £r(l)+ , £r(l)- 



For the process — > //g we have = i?^ + ffqq , withF] 



r $ + (Q 2 ,t,u,9,<l>) = V+^j^ + (t, u )-L(£o + ^o + ^2)-7e_(t,«)-^. 

S I IVlrp QrYlVIrp 



H$-(Q\t,u,eA) = ^{^((.^-^.(i-^A-W-^uJ^^a}. (77) 
Here, 

7e±(t,u) = (g 2 -t) 2 ±(g 2 -u) 2 . (78) 



16 The term proportional to >4i in Eq. (|7^) differs from the analogous expression in Eq. (16) of Ref. |19|j . 
In addition, Eqs. (16)-(18) of Ref. p9| are missing colour factors. 
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Similarly, for qg -> llq' (H$ = if «+ + we have 

#«+(Q 2 ,t, M ,M) = -v+^^(n + (t,s)-n + (o,s))c + 

- ((Q*-t) 2 + K-(t ,u))^-A 1 ]\ , 

9t J J 

H$-(Cf,t,u,e,<l>) = -V-^S [ (^n + (t,s)-n + (0,s)-2^-t(Q 2 -s))A 3 

- ^{2s(Q 2 -s)+>Jl + (t,s)} A 4 \ . 



K+(t,-s) (A0 + A2) 
(79) 



(80) 



Note that the crossing relationships for all other two-to-three processes are given as 

t,u,e,<t>) = H$(Q 2 ,u,t,n-e,<j>) , 

#<;>(Q 2 ,f,u,M) = H$(Q 2 ,u,t,n-8,<j>) , 
Hg\Q 2 ,t,u,9,<l>) = H$(Q\t,u,ir-8,n-(l>) , 
H^(Q 2 ,t,u,9,<j>) = H^(Q 2 ,u,t,9,7T-<j>) • 



(81) 



The corresponding expressions for and look simpler when written in terms 



of z A and z B , 



1 Q 2 

Z { qq(Q 2 , z A , z B , 0) = r q q(Q 2 , z A , z B ) — (V^C + V~gA 3 ) , 



£«(Q 2 ,z A ,z B ,#) = r qg (Q 2 ,z A ,z B )%(V^C + V; g A 3 ) 
Here, the functions functions r are given as 



(82) 
(83) 



r 'qq(Q 1 Z Ai z b) 



C F<^(Z 2 A + zl)5(l - T]Z A - T)Z B + Z A Z B ) 

e(Q 2 - q 2 T ) 5(1 - z B ) - Q(Q 2 - q 2 ) 5(1 - z A ) (i±|) 

(Q 2 \ : 



2Q(Q Z - q z T ) 6(1 - z A )S(l - z B ) 



Mr. 



(84) 
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and 



r qg (Q 2 , z A , Zb) = T R {z A (r] - z B ) z 2 A zl + (1 - z A z B ) 2 5(1 - r]z A -rjz B + z A z B 
- Q{Q 2 -Qt) [z 2 B + {l-z B ) 2 ]5{l-z A )} . 



(85) 



When using crossing relationships analogous to Eq. (|8l|) for £W one should remember 
that t «-> u implies z A ^ z B . 

From Eqs. & is evident that H {1) is an integrable function of q\ at q 2 — > 0. 

We now show this is also true of Eqs. ( |82| ) and fl83|). Define 

1 - r)X B # 1 - r?x A 



and consider the integral J a h 

/"^ dz & d/Zj: 

J a 



T] - X A 



X A Z A Jx B Z 



/a (—)/&(—) r ab (Q 2 ,z A ,z E 



•z k > 



Inserting the explicit form of r q q from Eq. (84) and using the representation 

q 2 f x A dz\ 

In -g = ln[?7 2 - 1] = In [(77 - x A )(r) - x B )] - \ 



lx A 7] — Z A 



we find 



J, 



[ dz A ( dz B 8(l - 7]Z A - 7]Z B + Z A Z B ) 

JXA Jx B 









f dz A 


~ J 


x A 


X 


9qq{ z Ai 






+ J 


f dz A 




x A 






+ J 


f dz B 




xb 



+ <? M -(M) 



9qq(z A , 1) - 9qq(h 1) 
g q q(l,Z B ) -g q g(l,l) 

(rj - x A )(rj - x 



0(x* A - z A ) &(Q 2 - q 2 T ) 



(V ~ z a) (1 - z A ) 

e(x* B -z B ) &(Q 2 -q 2 T ) 



(V - z b) (1 - z B ) 



111 



(1 - x A ){\ - X E 



e(q 2 - Q 2 ) In 



(1 - X A )(1 - X B ) 
7] 2 — 1 



with 



9qq(z A} Z B ) = C F — -fqi — jf, 



Za Z-p 



X B 
Zb 



(86) 



(87) 



(89) 



(90) 



26 



In this form the vanishing of J qq in the limit 77 — ^ 1 is manifest. Similarly, using Eq. (|35|) 



for r qg we obtain 



Jqg — / dZ- 



lx B 

where z* A = (1 — r]z B ) /(rj — z B ) and 



1 

9 q9 (z* A , z B )Q(x* B - z B ) - Q(Q 2 - ql)g qg (l, z B ) , (91) 



g gg (z A , z B ) = T R zlzl + {1 ZAZ * )2 f q (^)f 9 (^) • (92) 

z B Z A Z B 

Clearly, in the limit 77 — > 1 we have that J qg — > 0. 
A. 3 "Plus" distributions 

The definition of single plus distribution used in Section [I] and in this appendix is the 
standard one, 

1 dz f(z) [g(z)\, = I' dz [f(z) - f(l)]g(z) - f(l) f dz g(z) . (93) 







In Section |3] we also used the double plus distribution defined as 

/ dz A f dz B f(z A ,z B ) I- ^- -) = 

Jx A Jx B \(1 - z A )(l - z B ) ) ++ 

f 1 dz A C dz J {ZA > Zb) - f{ZA > 1} ~ Zb) + / 1} 

Jx A Jx B (1 — Z A )(1 — Z B ) 

+ f 1 dz A f(z A , 1) ( -J— ) ln(l - x B ) + ( l dz B f(l, z B ) {—^—) hi(l - x A ) 

Jxa V J- — Z A J _|_ Jx B \ 1 — Z B J _|_ 

- /(l,l)ln(l-x A )ln(l-x B ) . (94) 

B Numerical evaluation of Bessel transform 

When q and b become large it is convenient to use the asymptotic expansion^ to evaluate 
that portion of the Bessel transform because of the large cancellations between different 



17 This method was suggested in |l5j . The formula corresponding to Eq. (|9^) in Ref . |L5| contains an 
overall sign error. 



27 



cycles of the Bessel function. For large b we find that 

i 

dy Jo(qy)Ky) = 



x 



cos oo + n — I — — — nr 



Kb) 
b k+ l 



If we further make the special choice of 6 = b s 



2s 



the even n terms will cancel because 



1\tt 



4' q 



cos 



and we have 



dy Jo{qy)h{y) 



/ 7T 7T 

(go, + n- + - 



1 2 1 

9 V 7T(o6 s 



s = 1,2,3 



sin n— 

V 9 



b s h'{b s )--h{b s ) +0 



(95) 



(96) 



(97) 



(98) 



In practice it is found that the explicit integration needs to be done over only a few cycles, 
before approximating with the asymptotic expansion, Eq. (p8|). 
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Table 1: The range of values for g for which results shown in Figures |] through 10 were 
obtained. 



Experiment 


Range in Q [ GeV] 


Range in yjr 


Range in g [ GeV 2 ] 


R209 


5-8 


0.081 - 0.129 


0.52 - 0.55 


E288 


5- 6 

6- 7 

7- 8 

8- 9 


0.182 - 0.219 
0.219 - 0.255 
0.255 - 0.292 
0.292 - 0.328 


0.30 - 0.35 
0.32 - 0.40 
0.33 - 0.46 
0.38 - 0.50 


E605 


7- 8 

8- 9 
10.5 - 11.5 
11.5 - 13.5 


0.180 - 0.206 
0.206 - 0.232 
0.271 - 0.296 
0.296 - 0.348 


0.34 - 0.41 
0.42 - 0.49 
0.46 - 0.53 
0.46 - 0.51 



Table 2: Vector boson couplings in the notation given in the Appendix |A]. Lepton and 
quark charges are Q e = —1, Q u = 2/3, and Qd = —1/3. 



Boson 


k 


/r 




9k 


W 


9w/V2 





9w/V2 





Z 


9z(T 3 - Q e x w ) 




9z(T 3 - QfXyr) 
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Figure 1: Scales ji(b) and \(b) compared to b /b for Q = 5, 100 GeV. 
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Figure 2: The Sudakov form factor for Q = 5, 10 and 100 GeV. 
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Figure 3: S da 2 / dy/rdy distribution from E288 compared to theory multiplied by K 
0.83. 
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Figure 4: E da 3 /d 3 p distribution from E288 compared to theory with K = 0.75 (full line) 
and K = 0.83 (dashed line). 
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Figure 5: S da 2 / 1 d^frdy distribution from E605 compared to theory multiplied by K 
0.88. 
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Figure 6: E da 3 /d 3 p distribution from E605 (mass bins with Q < 9 GeV) compared to 
theory with K = 0.88. 
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Figure 7: -E 1 da 3 /d 3 p distribution from E605 (mass bins with Q > 10.5 GeV) compared 
to theory with K = 0.88. 
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Figure 8: Best x 2 /dof obtained by varying g for R209 data (5 GeV < Q < 8 GeV and 
q T < 3 GeV). For comparison, with the LY non-perturbative function we obtained x 2, /dof 
of about 1.3. 
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Figure 9: Best x 2 /dof obtained by varying g for E288 data. We used four data sets 
with q T < 2 GeV and Q below 9 GeV: 5 GeV < Q < 6 GeV, 6 GeV < Q < 7 GeV, 
7 GeV < Q < 8 GeV, and 8 GeV < Q < 9 GeV. Theoretical results were multiplied by 
K = 0.83. For these data sets the LY form of the non-perturbative function yields x 2 /dof 
of 6.5, 14.1, 22.9 and 27.5, respectively. 
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Figure 10: Best \ 2 /dof obtained by varying g for E605 data. We used data sets for 
7 GeV < Q < 8 GeV, 8 GeV < Q < 9 GeV, 10.5 GeV < Q < 11.5 GeV and 11.5 GeV < 
Q < 13.5 GeV (with q T < 2 GeV and K = 0.88). For these data sets the LY form of the 
non-perturbative function yields x 2 /dof of 2.2, 3.5, 16.6 and 17.3, respectively. 
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Figure 11: Comparison of the theoretical q T distributions for W = W + + W~ production 
at Tevatron (s/S = 1.8 TeV) with (D(a>s) perturbative calculation. These results were 
obtained with gaussian (&n m = 0.5 GeV^ 1 and g = 0.33 GeV 2 ) and LY form of F NP . We 
assumed BR(W -> ev) = 0.111. 



43 




Figure 12: Comparison of the theoretical (resummed plus finite) q T distributions for E288 
experiment (y/S = 27.4 GeV, 5 GeV < Q < 6 GeV and -0.27 < y < 0.33) with 0(a s ) 
perturbative calculation. These results were obtained with gaussian (by im = 0.5 GeV -1 
and g = 0.33 GeV 2 ) and LY form of F NP . 
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Figure 13: x 2 /d°f obtained by varying g for CDF W + + W~ (full line) and Z da/dq T 
data (dashed line) below q T = 45 GeV. LY form of the non-perturbative function yields 
X 2 /dof of about 0.6 and 0.5 for W + + W~ and Z data, respectively. 
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Figure 14: Theoretical prediction for CDF W + + W da/dq T data. These results are ob- 
tained using several different parton distribution functions, and with an effective gaussian 
form of F NP (g = 3.0 GeV 2 ,6 lim = 0.5 GeV" 1 ). 
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Figure 15: Theoretical prediction for CDF Z da/dq T . These results are obtained using 
several different parton distribution functions and with an effective gaussian form of F NP 
(g = 3.0 GeV 2 ,6 lim = 0.5 GeV" 1 ). 
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Figure 16: Transverse mass distribution for W + + W~ production at Tevatron. We used 
LY form (with their central values for the non-perturbative parameters), and also an 
effective gaussian form of F NP (with g = 3.0 GeV 2 and fe lim = 0.5 GeV -1 ). 
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Figure 17: p l T distribution from W production at Tevatron. We used LY form (with 
their central values for the non-perturbative parameters), and also an effective gaussian 
form of F NP (with g = 3.0 GeV 2 and b hm = 0.5 GeV" 1 ). 
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